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SUMARY 


The general equation^ governing the fixed- trim water landing of a 
straight-keel seaplane with a hull of arbitrary constant cross section 
are presented in such a form that the landing motions and loads are 
expressed in tenas of the steady-planing characteristics of the sea- 
plane. In order to verify the general validity of these equations, 
solutions are made for the water landing of a rectangular flat plate 
and are compared with experimental impact data. Calculated and experi- 
mental time histories of draft, velocity, and load are in good agreement. 
A survey is made of the available information on seaplane planing char- 
acteristics which is suitable for use with the analysis of the paper. 

• INTRODUCTION 


The National Advisory Committee for Aeronautics has undertaken an 
extensive program of theoretical and experimental research on hydro- 
dynamic impact loads in order to establish a more rational foundation 
for water-loading requirements for the design of seaplanes. Most of the 
results of this program to date are contained in references 1 to 6 . The 
development of the theory in these various papers usually proceeds sub- 
stantially as follows: First, a theoretical or semiemplrical analysis 

is made for the hydrodynamic forces acting during the two-dimensional 
impact of a body on a smooth water surface; the three-dimensional Impact 
and planing case is then treated by assuming that the fluid flow occurs 
primarily in two-dimensional planes oriented normal to the keel and 
applying an approximate over-all correction to account for the differ- 
ence between the two-dimensional and three-dimensional . cases . This type 
of approach to the impact and planing problems has been found to provide 
fairly reasonable estimates of the impact loads on certain types of sea- 
plane hulls (refs. 1 to 6), particularly those with scalloped bottoms 
and V-bottoms when the chines are not immersed below the water surface. 
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However, for impacts involving chine Immersion, including impacts of a 
rectangular flat plate, accurate two-dimensional solutions have not yet 
been derived for all cases. For such cases, or for cases where greater 
accuracy is desired than can be obtained from the two-dimensional analogy, 
other procedures for computing impact loads and motions must be developed. 
The purpose of this paper is to develop such’ a procedure by relating the 
basic seaplane Impact equations to the planing characteristics of a sea- 
plane and to present the solution of these equations in such a form that 
the impact loads and motions may be calculated from these planing char- 
acteristics. This derivation begins with the same assumptions as were 
made in the preceding papers on impact theory (refs. 1 to 6), namely, 
that the instantaneoiis forces dizring a landing and in planing depend 
only on the components of motion normal to the keel. The differential 
equation of motion based on these assimptlons is presented and then inte- 
grated to obtain equations for the time histo'ries of the draft, velocity, 
and hydrodynamic load during an Impact. The validity of the resulting 
equations is then tested by comparisons of experimental and calculated 
Impact loads and motions for water landings of a rectangular flat plate. 
Finally, a survey and evaluation is made of the available Information 
on the planing characteristics of seaplanes. 


SYMBOIS 


A hydrodynamic aspect ratio, (Wetted length at teel)^ 

Wetted area projected normal to keel 

t 

b beam of model, ft 

Cg function of angle of dead ris'e (O.75O for 3 = 30°) 

Ff hydrodynamic force normal to keel (normal to surface for a flat 

plate ) , lb 

F^ vertical hydrodynamic force, lb 

g acceleration due to gravity, 32.2 ft/sec^ 

H perpendicular distance between keel and plane of chines, ft 

m mass of model, slugs 

virtual mass, slugs 

n^^ Impact load factor, Fg/mg 



wave-rise factor (see fig. 1(c)) 

time after water contact, sec 

instantaneous resultant velocity of model, fps 

instantaneous velocity of model parallel to undisturbed water 
surface, fps 

instantaneous draft of model normal to undisturbed water 
surface, ft 

instantaneous velocity of model normal to undisturbed water 
surface, fps 

instantaneous acceleration of model normal to undisturbed water 
surface, ft/sec^ 

angle of dead rise, deg 

Instantaneous fll^t-path angle relative to undisturbed water 
surface, tan”^ f-, deg 

X 

perpendicular distance between keel and point of initial water 
contact, ft 

Instantaneous velocity of model normal to keel (normal to model 
svirface for a flat plate), x sin t + z cos fps 

Instantaneoiis acceleration of model normal to keel (normal to 
model surface for a flat plate), ft/sec^ 

wetted length at chines, beams 

z 

length of model "below imdistxirhed water surface^ h sin ^ 
beams 

wetted length at keel or longitudinal center line of model, 
beams 


^e ^ch 

mean wetted length, ^ , beams 

instantaneous velocity of model parallel to its longitudinal 
center line, x cos t' - z sin t, fps 



k 


MCA TN 2814 


p 

mass density of water, 1.938 alugs/cu ft 


fC 

T 

trim, deg 



cp(A) 

aspect-ratlq' correction 



ilf(aj) 

psi- function, ^ + log^ oi - 1 



SulDscripts: 


— 

ch 

at chine Immersion 



f 

rectanigular flat plate 



max 

msiximum value 



o 

% 

at water contact 



pi 

steady-planing conditions 



r 

at rebound 


i 

Dimensionless vetriables; 



Ca 

beam- loading coefficient, m/pb^ 




lift coefficient for steady planing, Fg/^Vpi^b^ 



• - 

Cb 



03 - 

l + 2w 
m 



Cl 

impact lift coefficient, 



Cy 

speed coefficient for steady planing (Froude number). 

Vpi/V^ 


k 

oz/b 

generalized draft coefficient, — / 

2C^ sin^T coa^'T'^O 

Cb’4 


e 

Impact parameter, tan(7o '’r}/tan t 



K 

impact approach parameter', cos (7r, + t ) 

sin 7o 
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Dots axe used to indicate derivatives with respect to tlmej for 
example. 


z = z = ^ 

dt dt 

ANALYSIS 

Derivation of Equation of Motion 


In the first part of the analysis, the basic differential equation 
for the hydrodynamic loads and motions occurring during the oblique 
water impact of a straight-keel seaplane with a hull of arbitrary con- 
stant cross section is presented and then converted into such a form 
that the equation is expressed in terms of the conventional planing 
coefficient Cg. Throughout this analysis the seaplane is assumed to 

remain at a fixed trim and to have zero rpll and yaw. The wing lift 
force acting on the seaplane is assumed to be equal to the wei^t of 
the seaplane , so that the net force acting to accelerate the seaplane 
is the hydrodynamic force. 

During the landing of a long narrow body of constant cross section 
on a smooth water surface, the hydrodynamic force on the body is 

generally assumed to be composed of two terms, one proportional to, the 
square of the component of the model velocity normal to the keel ^ (see 
fig. 1) and one proportional to the normal deceleration of the body 
(see, for example, ref. l), or 

+ e’^’ (1) 


where D and E are coefficients that depend on the instantaneous 
gometrical conditions (body shape, trim, and draft). The last term of 
equation (l) can be considered to represent the inertia force corre- 
sponding to the acceleration of the virtual mass of water associated 
with the impacting body; the coefficient E is thus equal to the virtual 
mass of water associated with the body (E = i%). A procedure for com- 
puting this quantity is given in appendix A. The coefficient D can be 
interpreted in terms of the conventional planing coefficient Cg as 

follows: In steady-planing tests the deceleration is zero and the 

normal^ velocity is related to the planing velocity Vp 2 _ by the rela- 
tion ^ = Vp 2 _ sin T (see fig. 1(a)). Also, the normal hydrodynamic 
force (neglecting friction) is related to the, vertical force by the 
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relation = F^ cos t. Substituting these relations into equation (l) 
gives for planing 


Fg. = DVp-j^^ sin^T cos t 


Since the planing coefficient C-g is defined as 


r* _ • 

® " I pVpi^bS 


( 2 ) 


(3) 


a comparison of equations (2) and (3) shows that D is related to Cg 
by the equation 


D 


■u2« 

pb Cg 


2 sln^T cos T 


(4) 


5k 


By use of this relation (and E = m^), eqiiation (l) can be rewritten as 


^5 = 




2 sin T cos T 


t 


(5) 


Since the wing lift on the model has been assumed to be equal to the 
weight of the model, .the hydrodynamic force is the net force acting on 
the model, and according to Ifewton’s second law this force must be equal 
to the inertia reaction, or 


F 




2 


pb Cg 

~~2 

sin T cos T 


%£ 


( 6 ) 


Equation (6) is the general differential equation for the force and 
motions for a seaplane landing, expressed in terms of the conventional 
planing coefficient Cg. 
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Solution of Equation of Motion 

In order to solve the differential equation of motion it is con- 
venient first to rewrite equation (6) in the form 



pb^ 


B 


2 sin^T cos T^l + 


2 


Then, converting from nonnal force to vertical force by multiplying 
through by cos t and abbreviating 




1 + — 
m 


gives the relation 


;• pb^C- 

Fz = -m^ cos T = L I 

2 sin T 


(7) 


( 8 ) 


The quantity C-g’ for . a given hull shape, weight, and trim depends 
on only one variable, the draft. On the other hand, ^ and ^ depend 
on the two variables ^ and tj equation (8) is thus an equation in 
three unknowns. In order to solve this equation one of these variables 
must be expressed in terms of the other two. The variable which is 
simplest to eliminate is which can be eliminated by the following 
substitutions (see fig. 1(b)): 

= — - — + \ tan' T ■ (9) 

cos T 


•• -7 

5 ^ 

COS T 

» 

where | is the component of the seaplane velocity parallel to the keel, 
which, from the previous assumptions of no friction drag and equal wing 
lift and float weight, is assumed to remain constant throu^out any 
landing. Eqiiation (8) can then be expressed in the form 

pb^^Cg 
2m sin t 


lCOS t 


+ I tan T 
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(z + i sin t) 


pb^B* 

2 = 2 2 

2m sin r cos t 


Since z dz = 
in the form 


— dz = ^ d& = z dZj equation ^(ll) can he also expressed 
dt dt 


z dz ph^Cg'dz 

2 2 2 
(z + i sin t) ' 2m sin t cos t 


z ^ z 

T i Bln T 

2 2 2 
' ' 2m sin t cos t 


I sin T 


which can he easily integrated to give 


IsinT T Zo + |sinT isinr z + lsinr 

^ + loge— ; T : l°Se “ : 

Zq + I sjLu t gsinr z+l sin t | sin t 


2 2 / 
2m sin T cos t ^0 


Cg dz 


The quantity on the ri^t-hand side of equation (l2), which depends only 
on the draft (for a given hull shape, weight, and trim), will henceforth 
he abbreviated as 


2 2 
2m sin T cos 


.r C3-. 

T^O 
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and will "be referred to as the generalized draft. An alternate form 
of equation ( 13 ), in terms of the heam- loading coefficient, is 


k = 


2C/ 


p 

sin T 


COS^T 




(li^) 


where 


Ca = 


m 


( 15 ) 


ph- 


In the remainder of this paper Cg is assumed to he a known function 
of draft which has heen determined either experimentally or theoreti- 
cally. Since all the other factors in equation (l4) are known con- 
stants, except the quantity m^/m which was used for converting Cg 

to Cg’ and which can also he expressed as a function of draft as indi- 
cated in appendix A, the quantity k can also he considered as a known 
function of draft. 

In order to compute the vertical hydrodynamic force from equa- 
tion (8), equation (12) must he expressed in terms of ^ according to 
equation (9) as follows: 


^o 


loSe r 




tan T 


- 1 - 


I tan T 




+ logg 


tan T 




I tan T 


- 1 = k ( 16 ) 


(The quantity 1 has heen arbitrarily added and subtracted from the left- 
hand side of this equation for later convenience. ) In order to facili- 
tate computation, a function i)r is defined as follows: 


= ^ + logg CO - 1 


( 17 ) 


Equation (l6) can he expressed in terms of the ilc -function as (compare 
eqs. (l6) and (l7)) 
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or 


tT- 


t 


tan Tj 






tan Ti 


(18b) 


or, in terms of the inverse function of as 




I tan T 


= 


-1 




^o 


tan T> 


(l8c) 


or 




-1 






tan TJ 


- k 


t 


I tan T 


(I8d) 


The constant Iq/1 tan t will he henceforth ahhreviatexi as g. Since 
I (which equals Iq) is equal to cot(7o ^ ^ 

can also he expressed in the form 




e = 


tan(7o + t) 


& tan T 


tan T 


( 19 ) 


It is noted that this constant € is related to the impact approach 
parameter k of reference 3 "by the relation e = 1 + i, where 


K 


sin T 
sin 7 q 


cos(7o 


+ t) 


The values of the i|f-function necessary for the use of equation (l8d) 
are given in table I and in figures 2(a) and 2(b). The procedure for 
using figures 2(a) and 2(b) is illustrated in figure 2(c). For example, 

to compute the quantity i|r"^[t(e) - first find e on the cs-axis. 

The corresponding ordinate is t(e). Then go- down on the Ir-axis a dis- 
tance k to obtain - k. The quantity ■'lf”^[^(6) - is then the 
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point on the co-axis corresponding to '>)f(e) - k on the ilf-axis. (For 
more accurate results a similar procedure can be followed with the aid 
of table I. ) It is noted that the inverse \|f-function is a double-valued 
function. The physical significance of this double-valued nature for 
the present problem is as follows: In each impact each value of draft 

(and of k) is reached twice, once while the seaplane is going down into 
the water and once while it is coming out of the water. Values of 



greater than unity correspond to the descending motion 


and those less than unity correspond to the ascending motion. 


The vertical hydrodynamic force is now obtained as a function of 
the generalized draft by combining equations (8), (l8d), and (I 9 ) to 
yield 


F 


z ~ 


,2 -.a, 2 ^ 

pb e tan t 

^ . 2 
2 sin T 




2 




( 20 ) 


The result may also be expressed as a dimensionless lift coefficient 



or as a load factor 




_ F^ _ P 


b\^cos.^(7o + t) 


mg 




2mg cos T 




( 22 ) 


The vertical velocity can be related to the generalized draft by 
combining equations ( 9 ), (l8c), and ( 19 ) to give 
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The rehouQd (and Initial) 
the generalized draft equal to 


It should he noted that the rehound velocity depends only on the impact 
parameter € and is not affected hy the shape or mass of the seaplane. 

The maximum draft is obtained hy setting the vertical velocity 
equal to zero in equation (23) so that >1^" ■^|t(e) - = 1 or 

■<|f(e) - k ^nny = ■'jf(l)^ and since ■>lr(l) = 0 (see table I), 

W = t(e) (25) 

where the relation between and z^„y is given hy equation (l4). 

The relation between the time and the draft is given hy 

t = r ^ dz 

Jq dz 

that is, a plot is made of the quantity l/z against z and is graphi- 
cally integrated to obtain the time. 

A proposed procedure for computing impact loads and motions 
according to the preceding equations is given in appendix B. 



vertical velocity is obtained by setting 
zero: 




e - 1 




Substantiation of Theory 

The results of the preceding derivation can be compared with experi- 
mental impact loads and motions in order to test the validity pf these 
results and thus of the assumptions upon which the derivation is based. 
For the non- chine -Immersed case the theory of the present paper is sub- 
stantially the same as the theory of references 1 to 6 ; hence, the agree- 
ment between theory and exj)erlment demonstrated in those papers for both 
the planing and the impact canes constitutes a verification of the theo- 
retical equations of this paper. However, for impacts which involve more 
than a slight degree'"of chine immersion, such as those of a rectangular 
flat plate, no substantiation of the theory has been previously given. 
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To supply such a suhstantiation, experimental arid calculated impact 
loads and motions for a rectangular flat plate are compared in figure 3. 
Experimental time histories of draft, vertical velocity, and hydro- 
dynamic load are shown for five landings. The planing data needed for 
these calculations were obtained from reference 7 (see appendix C for 
details of computations) and the corresponding impact data were obtained 
from references 8 and 9- Tiie good agreement shown in figure 3 for the 
calculated and experimental loads and motions substantiates the general 
validity of the impact equations. The small discrepancies that do exist' 
are clearly- within the limitations of accuracy of the experimental impact 
measurements. 


DISCUSSION OF PLANING DATA AND ANALYSES 


In view of the preceding substantiation of the impact' equations of 
this paper, the remaining problem is to provide reliable means of pre- 
dicting the planing characteristics of seaplanes. The available infor- 
mation on this subject can, for convenience, 'be grouped into the fol- 
lowing subdivisions: (l) the i^ctangular flat plate, (2) prismatic 

bodies with angles of dead rise greater than 10°, (3) prismatic bodies 
with angles of dead rise between 0° and 10°, and ( 4 ) prismatic bodies 
of arbitrary constant cross section. These various- subdivisions will 
now be considered in detail. 


The Rectangular Flat Plate 

A large quantity of flat-plate planing data are available in refer- 
ences 7 and 10 to 13. Hpwever, much of this information is Incomplete 
or was obtained at such low planing speeds that the buoyant forces (which 
are considered negligible during landings) are significant, and these 
data cannot always be safely extrapolated to find the high-speed 
(buoyancy- free ) planing characteristics needed for impact calculations. 
Also, these available data are usually limited to relatively low trims 
or small ratios of wetted length to beam. Consequently, it is doubtful 
that the available data are adequate to cover all practical landing 
conditions. 

The available analyses of flat-plate planing data, notably refer- 
ences l4 to l6, are all primarily empirical. They are consequently use- 
ful for Interpolating the experimental data, but they cannot necessarily 
be extrapolated accurately beyond the range of the experimental data. 

One difficulty in using these empirical analyses is the fact that they 
are restricted to the determination of the relation between the planing 
coefficient C-g and the ratio of wetted length to beam, whereas for 

Impact calculations the relation between Cg and the draft must be 
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known. In the case of the flat plate the water rises up in front of 
the plate (see fig. 1 (d)) so that the wetted length materially 

exceeds the length of the model below the undisturbed water surface 

where Xj is related to the draft, by the equation X .5 = . 

^ ^ b sin T 

Planing data showing the relation between Xg and X^j are plotted in 

figure hf together with a simple empirical equation fitting the data. 
This relation appears to be substantially Independent of the trim for 
the range of data available, namely, for 'trims below l 8 °. For higher 
trims, according to impact data in reference 9 , the wave rise Xg - X^ 

increases with increase of trim. (The quantity Xe used in this paper 
roughly corresponds to the quantity Xp in reference 9.) 


Prismatic Bodies With Angles of Dead Rise Greater Than 10° 

The problem of the determination of the planing characteristics of “ 
a V-bottom hull where the chines do not penetrate the water surface has 
been treated rather extensively, both theoretically and experimentally, 
in references 1, 3^ 5^ and 6 for angles of dead rise between 10° and 4o°. 
However, when the chines do penetrate the water surface the analysis 
becomes more difficult. The first important problem is to determine 
when chine immersion occurs or, more specifically, the draft at which 
chine immersion occurs Figure l(c) shows that this occurs when 

the splashed-up water crosses the chines, or the draft is 


_ H cos T b tan ^ cos t 
Zch - R 


(27) 


(In fig. 1(c) water pile-up in front of the hull has been neglected, or 
^d ~ ^e* ^ wave-rise factor R has been evaluated in different ways 

by several investiga,tors. According to Wagner’ s two-dimensional deri- 
vation for small angles of dead rise (ref. 17), R = ^. This value of R 

is supported by the theoretical solutions of reference l 8 , which also show 
that R decreases with increase of dead-rise angle, and by the experi- 
mental data of reference 15 for angles of dead rise between 10 ° and 30 °. 
Some full-scale landing data from a flying boat with 20° angle of dead 
rise (ref. I 9 ) indicate a value of R = 1.4. The impact theory of refer- 
ence 3 uses a value of R = - .l^tan 3 and reference 6 uses a value 

of R = 1. It should be noted that the last two values of ' R do not 
directly correspond to the physical water splash-up but are, rather. 
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effective values of splash-up. The first value of R was proposed from 
considerations of the over-all load prohlem and the second value resulted 
from the semiempirical analysis of planing data given in reference 6. In 
view of these observations it appears that R is approximately equal 
to rt/2; however, further research on this subject would be of value. 

Various analyses have been made for the planing coefficients of the 
V-bottom hull with chines substantially immersed, notably those of refer- 
ences 6, 15^ and 20. Reference 6 presents a semiempirical treatment of 
the problem which also considers the effects of buoyancy on the planing 
loads. In order to bring this treatment into agreement with the extremely 
limited data available at that time (refs. 6, 11, and 13) it was necessary 
to assume that the wave-rise factor R was equal to 1. Recently, more 
V-bottom planing data have been published (ref. 15)^ and from a consider- 
ation of these data as well as from the consideration of some impact data 
(to be discussed subsequently) for a float with 30° angle of dead rise, 
the analysis of reference 6 still appears reasonable for its two end- 
point conditions, the non-chine-immersed case and the deeply chine- 
immersed case, but it is somewhat conservative in the intermediate range 
as a result of the assumption that the wave-rise factor is unity. In 
reference 15^ where the more recent data are used, an empirical analysis 
restricted to the case of the chine-immersed hull is given which corre- 
sponds to the more realistic splash-up factor of ir/2 and which is- in 
fair agreement with the available experimental data.- (A somewhat parallel 
but less detailed treatment is given in reference 20. ) However, since 
the available data are still not very extensive and do not adequately 
cover the practical limits of large ratios of wetted length to beam, it 
is not yet certain that the empirical formula of reference 15 can be 
safely extrapolated to large ratios of wetted length to beam. There- 
fore, a need apparently exists for further study of the planing charac- 
teristics of V-bottom hulls with angles of dead rise greater than 10° 
for the case of considerable chine immersion. 

As a test of the relative merits of the two planing equations of 
references 6 and 15^ calculations based on these two formulas were made 
for the landing loads on a prismatic V-bottom hull having an angle of 
dead rise of 30° and are compared with experimental impact data in fig- 
ure 5« (since the analysis of reference 15 is restricted to the chine- 
immersed case, the planing coefficients prior to chine immersion were in 
both cases determined from the equation of reference 6 as is discussed 
in appendix D. ) At the low flight-path angle (/q = 2.2°) of figure 5(a), 
which corresponds to a landing in which the chines become immersed, but 
not deeply, the analysis of reference 15 gives much better agreement with 
the experimental impact data both in the shape of the load time history 
and in the magnitudes. This better agreement is due to the fact that 
reference 15 makes use of a more realistic estimate of the draft at which 
chine immersion begins than does reference 6. From an examination of fig- 
ure 5(b), which gives computed and experimental data for a large flight- 
path angle =5*9°) which corresponds to substantial chine immersion. 
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both analyses give reasonable qualitative representations of the load 
time history, reference 6 being somewhat conservative and reference 15 
being somewhat unconservative. Additional calculated and experimental 
impact load data shown in figure 5(c) substantiate these observations. 
In figure 5(<i) are shown calculated and experimental maximum drafts. 
These indicate the same conclusions as the load comparison j the maximum 
draft is underestimated by r&ference 6 (which is the same as over- 
estimating the load) and reference 15 tends to overestimate the draft. 

To summarize the preceding discussion of the planing coefficients 
for angles of dead rise greater than 10°, the use of the equations of 
references 6 and 15 apparently results in a reasonable calculation of 
the landing loads, reference 6 giving somewhat conservative results aind 
reference 15 giving- somewhat unconaervative results. However, for 
greater accuracy further planing research seems to be needed. 


Prismatic Bodies With Angles of Dead Rise Between 0° and 10° 

As yet no. experimental data are available to test the existing anal- 
yses for angles of dead rise bet-vfeen 0° and 10°. The analysis of refer- 
ence 6 is known to be inadequate for the case of 0° angle of dead rise; 
therefore, it is as yet somewhat doubtful that this eqiiation can be used 
in this range (0° < 3 < 10°). The equation of reference 15# by the 
nature of its derivation, can be reasonably assmed to apply for dead- 
rise angles in this range. However, one difficulty arises in the use of 
this eqiiation in this range. The results of reference 15 are essentially 
expressed in such a form that the planing coefficient is related to the 
wetted length at the keel (\e in fig- l(<i)) rather than to the draft. 

For angles of dead rise greater than 10°, reference 15 states that the 
relation between the draft and the wetted length at the keel is approxi- 
mately equal to the sine of the trim { Xe = ^d = — ). However, for 

\ b sin T/ 

the case of the flat plate (0° angle of dead rise) this relation was 
previously shown to be considerably more complicated because of the wave 
rise in front of the plate (see fig, 1(d)), Thus, if reliable computa- 
tions are to be made for the landing loads for V-bottom surfaces with 
angles of dead rise between 0° and 10°, a need exists for an investiga- 
tion of the wave-rise effects in this range. 


Prismatic Bodies of Arbitrary Constant Cross Section 

For prismatic hulls of arbitrary constant cross section, little 
experimental Information is available on planing characteristics (ref. 11). 
However, reference 2 presents a theoretical means of analysis which is 
applicable to hulls of constant cross section without reflex chines or 
chine strips, and which has been partly substantiated by some experimental 
Impact data in the same reference. 
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CONCLUDING REMABES 


The equations for the landing impact of a seaplane with a hull of 
arbitrary constant cross section have been presented in such a form 
that the landing loads and motions are expressed in terms of the planing 
properties of the seaplane. Comparisons of calculated and experimental 
loads and motions for the landing of a rectangular flat plate have shown 
that this analysis is valid for a greater range of conditions than was 
covered by previous papers. A survey and evaluation of the available 
information on seaplane planing characteristics which is suitable for 
use with the analysis of this paper has indicated a need for further 
research on these characteristics. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , July 7, 1952 
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APPENDIX A 
VIRTUAL MASS 


In reference 2 an approximate expression is given for the virtual 
mass of a body of arbitrary constant cross section prior to chine immer- 
sion. This equation, converted -to the terminology of the present paper, 
is 


m 


2C^ sin' 


-I p 

\ — r / '=B If 

> COS^T^O 


(Al) 


This relation, substituted into equation (l4), gives 

s/b i n_ 

/ 

k = 


2 — 
2C^ sin T cos T 


1 + 


1 r'^ 

2 t cos^t Jn 


2C^ sin^ cos 




(A2) 


prior to chine immersion. 

Subsequent to chine immersion the virtual mass m., associated 

^ch 

with the wetted area forward of the intersection of the chines and the 
water surface can be computed from equation (Al). The virtual mass of 
the region aft of this intersection can be approximated by a semicylinder 
of water circumscribed about the beam and corrected for the finite aspect 
ratio. The resulting equation for the virtual mass is 


"V = “Wch + ^ ^ch'P(A) 


(z ^ z„h) (A3) 
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% 

or 


m 




ch . pJtV 

8m 


m 


Xch9(A) = 


Ob ^ ef: 


2C^ sin^T cob^t'-'O 


(z 


m 


where cp(A) Is the aspect-ratio correction. For cases where the water 
does not splash up in front of the model, the wetted chine length 


(see fig. 1(d)) is equal to 


^ ~ ^ch 
b sin T 


and equation (Ai<-) becomes 


”^w 

m 


2C^ sin^T 


2 

cos T 




Jtq>(A) /z Zch ] 

8C^ Bin T\t 13/ 


= ^ch) 


(A5) 




For the flat plate the wetted chine length is the total wetted length 
= ^e) equation (Alj-) becomes 

^ = 8^ ^e»(A) ■ (A6) 

For arbitrary shapes with splash-up in front, the Tise of equation (A5) 
is suggested. The error introduced by neglecting the splash-up for this 
purpose will usually be small or negligible. 

According to Pabst (refs. 21 and 22) the virtual mass of a rectan- 
gular flat plate of aspect ratio A, as determined by vibration tests, 
is equal to that of a circumscribed cylinder times the quantity cp(A), 
where 


cp(A) 



(0 < A < oo) 


(AT) 


which can be approximated as 


9(A) 


1 - 


2A 


(A > 1.5) 
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These relations are suggested for use in the computations of the aspect- 
ratio correction for models of arbitrary shape, provided that the aspect 
ratio of the modeL-can he computed as a fimction of draft without too 
much difficulty. If this procedure is difficult, however, the aspect- 
ratio correction may he considered approximately equal to unity: 

cp(A) % 1 (A9) . 

For cases where A is reasonably large, no significant error should 
result from the substitution of equation (A9) into equation (A5). 

The hydrodynamic aspect ratio of a model is, by definition, equal 
to the quotient of the sqiiare of the keel wetted length and the pro- 
jected wetted area normal to the keel. Thus, for a rectangular flat 
plate, 

A = Xq (AlO) 

and, for a wedge of dead-rise angle 3 with no forward splash-up of 
water. 



z tan 3 COB t 


b to 


% 


(All) 



t 
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APPENDIX B 


PROPOSED PROCEDURE FOR COMPUTATIONS OF IMPACT LOADS -AND MOTIONS 


In order to utilize the equations of this paper to compute the 
water loads and motions of a seaplane, the following procedure is 
suggested: 

(lO Select values of C^, t, 7^, and V^. 

(2) Make a plot of C-g against z/b. This plot may be obtained 

from experimental planing data or by utilizing theoretical planing 
equations. 

(3) Graphically or analytically integrate the Cg curve to obtain 
>z/b 

the quantity / C-r as a function of draft. 


fo 


(l^.) Compute nv/m as a function of draft according to equations (Al) 
and (A5). (For sll^tly better accuracy in the case of the rectangular 
flat plate, use equation (a 6), together with the empirical ciurve.in fig- 
ure If-, instead. ) 

(5) Compute Cg’ as a function of draft according to equation (7). 

(6) Compute k as a function of draft according to equation (l4). 

(7) Compute e from equation (I9) and then i|r(6) from table I 
or figure -2. 

(8) Select a series of values of k between 0 and ^(e). 

(9) For each value of k (the generalized draft), compute 

^ ^ as indicated in figure 2(c), with the aid of table I or 

figure 2(a) or 2(b), ^The larger value of t"^[t(6) - ^ obtained 

corresponds to the seaplane descending into the water and the smaller 
value corresponds to the seaplane rising out of the water. ^ 

(10) Compute the vertical force from equation (20)., (21), or (22). 
Plot this quantity against draft. _ The maximum load can then be obtained 
from this curve. 
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(11) Obtain the vertical velocity as a function of draft from 
equation (23). 

(12) . Obtain the relation between draft and time by graphical inte- 
gration of equation (26). 
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APPENDIX C 

DETAILS OF CCMPUTATIONS FOR FIAT PLATE 


The impact loads and motions shown in figure 3 for the rectangular 
flat plate were computed according to the procedure of appendix B. The 
relations between Cg and z/h needed for step (2) of that procedure 

were obtained as follows. The relation between the wetted length Xg 
and the planing coefficient was obtained from a cross plot of 

Locke's high-speed planing data from reference 7 (see fig. 6), The 
relation between the draft and the wetted length was obtained from fig- 
ure 4, which is a plot of data obtained by Sottorf (ref. 10) and 
Sambraus (ref. 12). The empirical formula shown in this figure was used 
for all computations. The relation between and z/b was obtained 

by combining figure 6 with the empirical formula from figure 4. 
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APPENDIX D 

DETAILS OF CCMPUTATIONS FOR 30° ANGLE OF DEAD RISE 


The Impact loads and motions shown in figure 5 for a model with 
30° angle of dead rise were computed according to the procedure of 
appendix B. The relation "between Cg and z/b needed for step (2) of 

that procedure was obtained from the following equations taken from 
references 6 and 15: 

According to reference 6 (if buoyancy be neglected), 


Cb = 1.1*2* cot% sin t(i - (* < ^ch) 


Cg = 0. 355 jt sin t cos^t 




Kt) 

n 4 rr-^z tan 3 cos t\ 

Cg sin T cos J 


> (Dl) 


(Z > 2ch) 


R = 1 


where Cg is a function of dead-rise angle tabulated in reference 6 
(0.750 for 3 = 30°). 


According to reference 15, for the chine-inmiersed case only (if 
buoyancy be neglected), 



R 



0.006530^ 


0, 6 




(2 >^ch) 


(D2) 
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where 


Cb^ = 0.012Uk^f^r^'^ 
= ^(^e + ^ch) 




® h sin T 


■\ _ ■), tan B 

ch ® ” It tan T 


Since equation (D2) is not valid for the case of no chine immer- 
sion, in order to use this equation for impact computations equation (Dl) 
was used to predict the pl aning coefficient prior to chine immersion. 
Chine immersion was taken to occur at the draft at which equations (Dl) 
and (D2) predicted the same planing coefficient. This procedure for 
determining chine immersion was chosen only because for the specific 
case considered (3 = 30°^ t = 15°) this value of draft corresponded to 
a wave-rise factor of roughly it/2, which is in agreement with the value 
predicted by equation (D2). 
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y(o)). 



(c) Illustrative procedure for computing the inverse ^-function. 


Figure 2.- Concluded. 
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Figure 4.- The increase in wetted length due to wave rise for a rec- 
tangular flat plate during steady planing. “ * 
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Experimental data obtained from tests desci^ibed in reference 23 
la^ct ccKpatations based on analTsis of reference 6 
In^ot coBQ}utations based on analysis of reference 15 


Time, t, see 

I Impact load tiane biatory. 7q = 2.2P} 

Zq = 3*^ feet per second; = 89.3 feet 
per second. 


Tine, t, see 

(b) Impact load time history. 7 ^ = 5 - 9 °; 

Zq = 9 »*t- feet per second; Xq = 90. 9 feet 
per second. 
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Flight-^ath angle, Yq> 


(c) ifeiximum loads. 


(d) Maximvm drafts. 


Figlire 5*- Comparison of calculated and exjjerimental loads and drafts 
for a prismatic model with an angle of dead rise of 30°, a heam- 
loading coefficient of I 8 . 8 , and a trim of 15°. 
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Figure 6.- High-speed steady-planing data. 3 = 0°. 
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